Kink Dynamics in a Topological 4 Lattice 
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It was recently proposed a novel discretization for nonlinear Klein-Gordon field theories in which 
the resulting lattice preserves the topological (Bogomol'nyi) lower bound on the kink energy and, 
as a consequence, has no Peierls-Nabarro barrier even for large spatial discretizations (h ~ 1.0). It 
was then suggested that these "topological discrete systems" are a natural choice for the numerical 
study of continuum kink dynamics. Giving particular emphasis to the </> 4 theory, we numerically 
investigate kink-antikink scattering and breather formation in these topological lattices. Our results 
indicate that, even though these systems are quite accurate for studying free kinks in coarse lattices, 
for legitimate dynamical kink problems the accuracy is rather restricted to fine lattices (h m 0.1). 
We suggest that this fact is related to the breaking of the Bogomol'nyi bound during the kink- 
antikink interaction, where the field profile loses its static property as required by the Bogomol'nyi 
argument. We conclude, therefore, that these lattices are not suitable for the study of more general 
kink dynamics, since a standard discretization is simpler and has effectively the same accuracy for 
such resolutions. 

PACS numbers: 05.45.Yv; 47.11. +j; 05.10.-a; 11.27. +d 



Coherent structures in nonlinear Klein-Gordon field 
theories are an active topic of research in many branches 
of physics, arising cither as a fundamental excitation in 
particle physics or as a collective excitation in effective 
field theories. Of particular interest are the so-called 
topological solutions, where the "kinks" are probably 
among the most studied ones (see, e.g. |Q and Q). 

Focusing on the numerical study of kink dynamics, 
Speight and Ward || suggested a new discrete sine- 
Gordon system which preserves the so-called "Bogo- 
mol'nyi bound" on the kink energy 0], being also ap- 
plied in a (ffi system H and recently generalized to non- 
linear Klein-Gordon equations in one spatial dimension 
H . These systems were then called "topological discrete 
systems" (TDS) by the authors (since the Bogomol'nyi 
bound depends on the topological charge, i.e., on the 
boundary conditions at infinity), and they showed that 
such lattices do not present the so-called Peierls-Nabarro 
barrier, which, for continuum field theories, is a rather 
artificial and undesirable effect induced by standard dis- 
cretizations. 

By means of numerical experiments with the velocity 
of a single kink, these authors also showed that the nu- 
merical drift in a TDS is significantly reduced even for 
large spatial discretizations h w 1.0, and concluded that 
this construction provides an efficient and natural choice 
for numerically simulating kink dynamics in nonlinear 
Klein-Gordon models. If the latter statement were in- 
deed true, it is clear that we would benefit a lot in both 
computational time and memory usage by using coarse 
lattices with h « 1.0 when studying complex kink behav- 



ior in these field theories, such as breather formation and 
multiple kink (antikink) scattering. 

However, and for our disappointment, we found that a 
TDS choice for such dynamical problems is good only for 
fine spatial discretizations h sa 0.1, where a simpler stan- 
dard discretization is also accurate. Our proposal in this 
report is to show how we have investigated these facts. 

The equation of motion for our problem is the well- 
known Klein-Gordon equation: 



d tt (f> = d xx <f>-V' (<f>), 



(1) 



with subscripts denoting partial derivatives and primes 
denoting derivatives with respect to <p. The potential we 
have adopted is the (double-well) 4> 4 one: 



V(0 = -(1-^) 



(2) 



Notice that it differs from that of Speight || by a mere 
factor of 2. We have, of course, corrected the TDS equa- 
tions in order to account for this fact, and a test of va- 
lidity of this rather simple correction will be given below 
for the case of a free kink (All physical quantities in this 
work will be given in dimcnsionless units. For specific 
applications, a simple rescaling of the variables in the 
above equations can introduce the desired units). 

In order to integrate Eq. (|j) we use the standard stag- 
gered leapfrog scheme, which ensures second-order accu- 
racy in time (see, e.g. M): 
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(3) 



This should be contrasted to the simpler expression from 
a standard discretization: 



where superscripts (subscripts) denote temporal (spatial) 
indices, overdots represent partial time derivatives and 
primes are derivatives with respect to <p. The spatial 
derivative is evaluated with a standard second-order cen- 
tered difference [d xx (f>i ~ h~ 2 {4>i + i — 2<f>i + (f>i-i)] which 
coincides with the resulting discretization required by our 
particular TDS choice || . Being a symplectic algorithm, 
the leapfrog scheme is suitable to Hamiltonian systems 
(see e.g. ||) and therefore we expect good energy conser- 
vation throughout our simulations. Indeed, in a typical 
simulation with h = 0.1, the energy is conserved to better 
than one part in 10 3 . 



V'{ ( jy l ) = 4>i-^ 



(6) 



However, the main feature of the TDS potential is that 
it makes it possible to construct directly from the lattice 
theory the Bogomol'nyi lower bound on the kink energy: 



< E P 



h 

71 



(7) 



with A the forward difference operator Afi = (fi+i 
fi)/h and Ep the lattice potential energy given by: 
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FIG. 1. Time evolution of the free kink velocity for 
h = 1.2. The upper curve is for the TDS lattice whereas the 
lower one is for a standard discretization. The initial velocity 
of the Lorentz boosted kink is v = 0.6. 

The "standard" discretization of the potential (Q) is 
straightforward, and involves simply the replacement of 
the continuum arguments x and t by their discrete coun- 
terparts, e.g. i and n. Put briefly, the TDS approach 
adopts a rather different lattice potential, which for the 
</> 4 case was suggested to have the form (5) : 



VTDs(0i 



(4) 



with the temporal indices equal to n. Of course, in the 
limit h — > both lattice potentials are equivalent to the 
continuum one (||). Notice that the resulting equation 
of motion involves a field derivative of VrDs(0j) which is 
somewhat intricate ||: 



E f 



i(A^) 2 + F T Ds(^) 



(8) 



The kinks are static solutions of the field equations with 
energy equal to the above lower limit, namely 2\/2/3. 

The first obvious thing to be done is to reproduce the 
main numerical results of [^) in order to check our TDS 
implementation. This is shown in Fig. |l| for the case of a 
spatial discretization h = 1.2 and initial velocity v = 0.6. 
This result is remarkable, and it is tempting to use this 
lattice in other dynamical kink problems, as suggested in 
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FIG. 2. Kink position after a collision with an antikink 
(the collision occurred approximately at t ~ 23). All the three 
curves were obtained from a TDS lattice. The most accurate 
one is, of course, for h = 0.1, where a standard discretization 
gives the same plot. The final kink velocity, as obtained by 
taking the slope of the curves for t > 100, are 0.2777, 0.3049 
and 0.5654 for h = 0.1, 0.4 and 0.8, respectively. 
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Since we have restricted our study to kink-antikink 
problems (including bound states, such as breathers), we 
benefit from the symmetry of system [<p(x, t) — 4>(—x, t)] 
so that we need to integrate only half the original lattice. 
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We have adopted the "adiabatic damping method" of 
Gleiser and Sornborger || for the study of kink-antikink 
bound states, and a moving boundary with constant ve- 
locity equal to that of light (i.e., with v = 1 in our dimcn- 
sionless units) for the study of kink-antikink scattering. 
These two methods are very effective from the computa- 
tional point of view, in that they avoid the reflection of 
practically all the outgoing radiation by the boundaries 
without the need of using huge lattices. 
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FIG. 3. Relative energy error for the adiabatic damping 
method of Ref. || in comparison to an outgoing boundary 
for our breather problem. The data for t < 100 is not shown 
since AE = for this case. The parameters, in the notation 
of the above reference, are k = 0.005, po = 150 and L = 250, 
where 2L is the physical size of the lattice. 

By using as a guideline the very precise work of 
Campbell et al. ||, we have investigated the behav- 
ior of kink-antikink scattering above the critical velocity 
v c = 0.2598..., where there are no resulting bound states. 
In Fig. ^ we show a particular sample of our simulations, 
for the case of an initial kink velocity v = 0.4 (similar 
results are obtained for different v). The final kink ve- 
locity as calculated in the above mentioned work for this 
same initial velocity is Vf ~ 0.28, where the authors have 
adopted a very fine lattice, namely h = 0.01. Compar- 
ing this result with ours, we notice that our simulations 
are quite accurate even for h = 0.1 (see notes in Fig. |^). 
From this figure we also clearly see that the TDS lattice is 
accurate only when h w 0.1, where a standard discretiza- 
tion also works very well (it is not shown in this figure 
since it coincides with the TDS curve). Therefore, even 
though the TDS lattice has achieved a very good constant 
kink velocity in coarse lattices, it overestimated the final 
kink velocity during the kink-antikink interaction (notice 
from Fig. that the kink velocity is practically constant 
for t > 23). We suspect that this fact is related to the 



breaking of the Bogomol'nyi bound, since, during this in- 
teraction, the field configuration loses its static property 
which is required by the Bogomol'nyi argument. 

Another application of interest within the context of 
kink dynamics is the formation of kink-antikink bound 
states, resulting in the so-called breathers. Perhaps the 
most popular breather is the one in the sine-Gordon the- 
ory ||], being mainly characterized by an exact periodic 
oscillation of the field such that the resulting coherent 
structure is stable, i.e., it retains the energy localization 
even when t — > oo. For the 4 theory in one spatial di- 
mension, however, it is today known that no such stable 
structure exists pp| . However, an extremely long-lived 
structure very similar to a legitimate breather exists, and 
has been a topic of active research not only in one [PTl, 
but also in two and three spatial dimensions |8|JlT 
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FIG. 4. Breather energy for three different lattice spacing 
ft in a TDS system. Again, the h — 0.1 case coincides with 
the result of a standard discretization. 

Using as a reference the work of Geicke on the loga- 
rithmic decay of 4 breathers with energy E < 1 ]ll| ] , we 
have also numerically investigated the effectiveness of the 
TDS approach and, as already anticipated above, the re- 
sults are disappointing. The breather, as obtained in the 
above work, arises from a dynamical kink-antikink inter- 
action, with the kinks initially at rest and at a mutual 
distance equal to 1.6. Again, and as a first step, we have 
checked our numerical routines by comparing our results 
with the ones available in the reference work. The main 
observable adopted here is the energy within a certain re- 
gion around x = 0, i.e., the effective breather energy. By 
observing the same region as Geicke, namely —b < x < b 
with b — 100, and integrating up to t = 2000, the re- 
sulting energy was E = 0.9732, whereas Geicke's was 
E = 0.9812. 1 We suspect that this difference is related 
to the presence of some dim radiation reflected by the 



1 Geicke reported two values of the energy: one corresponding to the configuration where (f>(0, t) assumed a maximum and the 
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boundaries of Geicke's lattice [which, although alleviated 
by his boundary conditions, is still evident from Table I 
in his work, for which a greater value of the lattice size a 
would probably follow the observed tendency to decrease 
E(t = 2000) for fixed b = 100]. As already stated previ- 
ously, we have adopted a very efficient damping method 
to avoid this sort of radiation, and have found a set of 
parameters for which the method gives a very accurate 
result (better than one part in 10 4 , see Fig. |^) in compar- 
ison to an outgoing boundary (which never reflects the 
radiation, see ||). 

In Fig. |]we present another sample of our simulations, 
now for the breather energy in a TDS lattice for three val- 
ues of h. We feel that this plot is sufficient to show that 
a topological lattice gives very different results even for 
relatively fine lattices, such as h = 0.4. For more coarse 
lattices, such as h = 0.8, the results start to differ not 
only in the magnitude of the energy, but also in its profile. 
Of course, for larger still spatial discretizations the error 
grows without bounds and therefore the TDS approach 
is definitely not suitable for such problems. Again, we 
believe that this error is due to the breaking of the Bogo- 



mol'nyi bound since, although the kink and the antikink 
are static solutions of the field equations separately, the 
same is far from being true for their bound state, i.e., for 
the breather case. 

We conclude this report by remarking that for all the 
problems tackled here where the TDS approach is ac- 
curate, i.e. where h rs 0.1, a standard discretization 
is not only simpler, but also computationally more ef- 
ficient since it involves fewer operations (compare Eqs. 
(||) and (g), for example). We cannot disagree with the 
fact that such topological lattices are indeed a remark- 
able discretization for the study of free kinks on a lattice. 
However, in the face of the numerical experiments pre- 
sented here, we can not agree that this construction is 
an efficient and natural choice for numerically simulat- 
ing kink dynamics in nonlinear continuum Klein-Gordon 
models |,§. 
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other when (j>(0,t) assumed a minimum. The difference, however, is quite small (~ 0.004), so here we quote only the smaller 
one. 
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